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Abstract 



A general stochastic model is developed for the total interference in wideband systems, denoted as the 
PNSC(a) Interference Model. It allows one to obtain, analytic representations in situations where (a) 
interferers are distributed according to either a homogeneous or an inhomogeneous in time or space 
Cox point process and (b) when the frequency bands occupied by each of the unknown number of 
interferers is also a random variable in the allowable bandwidth. The analytic representations obtained are 
generalizations of Cox processes to the family of sub-exponential models characterized by distributions 
from the a-stable family. We develop general parameteric density representations for the interference 
models via doubly stochastic Poisson mixture representations of Scaled Mixture of Normal's via the 
Normal-Stable variance mixture. To illustrate members of this class of interference model we also develop 
two special cases for a moderately impulsive interference (a = 3/2) and a highly impulsive interference 
£SJ , (a — 2/3) where closed form representations can be obtained either by the SMiN representation or via 

function expansions based on the Holtsmark distribution or Whittaker functions. To illustrate the paper 
we propose expressions for the Capacity of a BPSK system under a PNSC(a) interference, via analytic 
expressions for the Likelihood Ratio Test staistic. 

Keywords: Interference models; Cox Process; Doubly Stochastic Poisson Stable Process; Isotropic a-stable; 
Complex a-stable. 



I. INTRODUCTION 

Modeling interference in ad hoc or cognitive networks [1], [2], [3], [4], [5], [6] is an active 
research area in wireless communications. A key attribute often present is that interference in 
modern wireless networks exhibit an impulsive nature, evidenced by the growing literature on 
the understanding and study of such features in wireless communications. 

Several modeling approaches for impulsive noises have been proposed, especially for impulse 
radio ultra wide band communications (Laplace, Generalized Gaussian, Cauchy, a-stable, Mid- 
dleton class A, an overview is, for instance, given in [7] and the references therein). Each of these 



models shares something in common: they are members of what one can define as heavy tailed 
or impulsive noise processes; each is a member of the sub-exponential family of interference 
distributions. In this sub-class of interference models there is one particular sub-family, the in- 
stable distribution, which enjoys a very rich and extensive literature on its characteristics and 
attributes both theoretically as a generalization of the Central Limit Theorem and in applied 
settings. Several papers have dealt with such models, see for example in the signal processing 
literature the discussion in [4], [5]. 

A. Brief summary of the -paper. 

In this paper we extend the family of a-stable interference models, proving that one can obtain 
and characterize a total interference model in the class of doubly stochastic Poisson-Gamma- 
complex Isotropic Stable Cox processes, which we refer to as the PNSC(a) class of models. This 
involves three main contributions: 

• First we extend the model, for instance presented in [2], to a doubly stochastic Poisson 
random field. Practically, this means that not only the number of interferers is a random 
variable but also the number of sub-carriers occupied by each interferer from the total set of 
available carriers is random. There are several applied settings in which such an approach is 
required, for instance in cognitive networks where secondary users would use several sub- 
carriers to adjust their desired bit rate but not all of the possible sub-carriers available. In 
addition, in this context, a strong contribution of the paper is to also allow the Poisson field 
to be inhomogeneous in time or space. For example, we include illustrations of the cases 
where users are sparser when further from the access point or when directive antennas are 
used so that interferers are only located in a sector of the plane. 

• Having derived and interpreted the resulting class of impulsive total interference models, 
labeled the PNSC(a) family of interference, we need to develop tools for practitioners to 
utilize these models. Therefore, the second contribution involves developing analytic and 
closed form representations. We note that a challenge of stable distributions is that, in general, 
they do not admit analytic expressions for the density and distributions. This makes, for 
instance, receiver processing difficult to adapt to such models. We overcome this important 
challenge utilizing several key features of the family of complex isotropic stable distributions, 
which we extend to the doubly-stochastic Poisson-Gamma-complex Isotropic Stable Cox 
process setting derived. Namely, the key attributes we utilize are that stable random vectors 
are closed under convolution, they admit scaled mixture of Gaussian representations under 
projection and they admit infinite series expansion representations under projection. We 
apply these in a general context and demonstrate their attributes for interference modeling 
under truncation approximations of such series. 

Then for two special cases we develop some closed form solutions resulting from special 
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function representations, which do not require infinite series or scale mixture representa- 
tions. The first one considers situations with moderately impulsive noise (a = 3/2, which 
corresponds to a channel attenuation coefficient a = 2.66 if we refer to [2], [3] for instance) 
and the second to a highly impulsive noise (a = 2/3 corresponding to a = 6). We demonstrate 
and discuss how such analytical results can be very useful in the understanding and design 
of nodes or network in impulsive interference. 

• Thirdly we derive novel expressions for the Likelihood Ratio Test statistic and study the 
capacity of a binary input soft output symmetric memoryless channel for the class of PNSC(a) 
interference models. 

B. Organization and Contribution 

The paper is organized as follows: in section[[TJwe provide a brief technical review of theoretical 
properties of a-stable random variables and vectors that will be of direct consequence for the 
proofs and derivations of results in this paper. In section [HI] we describe the system and the 
different assumptions that we make to obtain the theoretical results. Finally in section [TV] we 
derive the interference distributions in the different cases and provide examples to illustrate the 
proposed model. 

Below, we highlight some specific contributions involving two generic frameworks which 
are developed based on analytic representations of distributions corresponding to the "total 
interference" across transmission bandwidth in a wide-band wireless communication system in 
which the number of users transmitting is treated stochastically and the bandwidth they occupy 
is also treated as stochastic. This involves the following key contributions: 

• In Theorem [TJ we extend the representation of [8], [2], [3], [4] to derive the log Character- 
istic Function (CF) for the total interference at a given frequency, for a random number of 
potential interferers in a spatial region of transmission. We prove that the resulting log CF 
can be represented by the family of isotropic bivariate a-stable distributions. In doing so, 
we generalize existing results in the following ways: (a) we allow the number of users to be 
distributed according to three possible scenarios involving homogeneous Poisson (Model 1), 
temporally inhomogeneous Poisson (Model 2) and spatially inhomogeneous Poisson (Model 
3) point processes; (b) in Lemma [IT] we detail analytic distributional representations for the 
bivariate isotropic stable distribution based on a Scaled Mixture of Normals (representation 
1) and an exact Projection based univariate representation via a Cramer-Wold decomposition 
(representation 2). 

• In Theorem |2] we utilize the isotropic a-stable representations of the total interference at a 
given frequency (Theorem [TJ), to extend these results to derive analytic solutions for both 
the distribution function and density function of the resulting inhomogeneous (spatially or 
temporally) stochastic Poisson compound processes that models the practically important 
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total interference across the entire transmission bandwidth. This involves development of 
analytic solutions for the distribution of a compound process (Poisson mixture) of bivariate 
stable components and utilization of the closure under convolution of the stable components, 
based on work of [9]. 

• In Theorem [4] we generalize the results derived in Theorem |2] to analytic representations of a 
Cox process (doubly stochastic compound Poisson-Stable mixture) for the total interference in 
the context in which the mean occupied bandwidth of each potential interferer is generalized 
to a stochastic model. This allows for potential scenarios such as different occupancy average 
requirements per user; or time varying occupancy requirements. 

• In Theorem [7] we derive expressions for the Capacity of the PNSC(a) interference models 
expressed analytically according to a likelihood ratio test statistic. We derive several analytic 
expressions for the LRT in the case of PNSC(a) interference models, making evaluation of the 
capacity highly efficient for generally intractable stable models. We illustrate their accuracy 
as a function of the stable tail index. 

C. Notation 

The following notation is used throughout: random variables are denoted by upper case letters 
and their realizations by lower case letters. In addition, bold will be used to denote a vector or 
matrix quantity, upper subscripts will refer to a specific interferer and lower subscripts to the 
element of a vector or matrix. 

II. Background on Statistical Interference Modeling 
A. Interference modeling 

Interference from undesired active users in a network will be a strong limitation in future 
networks performance. The interference model has been studied widely in information theory 
[10], [11], [12], [13]. If the exact capacity is not known some close approximations have been 
derived. The question on how to deal with interference is however still an open problem. In this 
regard a lot of work on multiuser detectors for instance have been proposed [14] but also, more 
recently, some new schemes for interference alignment [13] or amplifying interference [15] have 
been considered. However, those works aim at avoiding the interference and generally require 
some costly channel learning mechanisms or synchronization techniques. 

An alternative perspective is to consider that a certain amount of interference will be un- 
avoidable. Under such an assumption, a robust interference model can allow an effective design 
of receivers and networks to limit the resulting impact of such interference. This is a powerful 
tool to study for example outage probability or connectivity in networks. For instance several 
works on stochastic geometry are based on similar interference models as we present in this 
paper [16], [17]. If we consider k,r to be a random variable representing the number of active 
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interferers. In a rather general framework, the total interference is a random variable expressed 
according to Y = J2k=i^ipk where (ipk)k=i,-,N are independent, identically distributed and 
bounded random variables with even probability density function that depends on the physical 
layer design (see [1], [2], [3] for different examples). The (A k ) k=1: ... :N are positive, independent, 
identically distributed random variables that depend on the channel characteristics and determine 
the statistical properties of the total interference Y. 

To proceed, the most intuitive statistical approach would be to consider the asymptotic be- 
havior of the distribution of the total interference Y ~ F(y) and to determine under what 
conditions such an interference would belong to the domain of attraction of a Gaussian family 
of distributions, denoted T>q(F). Such an approach involves considering an asymptotic regime 
where the number of interferers grows to infinity while the contribution of each interferer to Y 
becomes infinitesimal. In non-impulsive, non-sub exponential distributional settings, this would 
typically result in application of a form of the celebrated Central Limit Theorem: Y converges 
in law to a Normal distribution, such that F(y) £ V G (F). However, in the general case in which 
impulsive noise is present it is well known that this asymptotic regime is not easily reached 
(see for example an in-depth study in [18] for impulse radio ultra wide band signals). Instead, 
the domain of attraction of impulsive noise models from the sub-exponential family, which are 
convolved to create the total interference Y ~ F(y) can belong to the domain of attraction of a 
stable family of distributions, denoted V S (F), for which V G (F) C T>s{F). 

A common requirement for convergence of such a sequence of i.i.d. interferences to converge 
to the Gaussian domain of attraction, involves a restriction on the variance of such summands in 
the sequence. This is not present in impulsive noise processes. One could argue that this feature 
may seem natural since it represents a channel attenuation, which by its very nature must be 
finite. However, the interference which is being modeled is actually compared to the desired 
link attenuation and can, in comparison, be "very large" and impulsive in nature. Such large 
impulsive realizations of the interference happen infrequently in practice but are sufficient to 
give an impulsive nature to interference. To capture these situations, heavy tailed distributions 
with infinite variance can be well suited while models with finite second order moments will fail 
to adequately capture such impulsive attributes observed. The generalized central limit theorem 
has then to be used (see [19, p. 22] or [20, p. 9]) and states that interference (for large kr) falls 
in the domain of attraction of a random variable with a stable distribution, F(y) £ V S (F). 

A general framework is proposed in [4] and application to cognitive radio with a modified law 
(truncated a-stable although the term truncated is slightly misleading) is presented in [6]. Here 
the truncation refers to a form of soft "tempering" of the stable distribution tails, as opposed to a 
hard thresholding. To prove the validity of the ct-stable assumption, the usual solution is to write 
its log CF as ifY (oj) = —a\tu\ a . This can be done in many situations (users' repartition, channel 
conditions, physical layer etc., refer to [2] for more details). One strong advantage of this model 
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over other proposed solutions is its theoretical foundations which we derive in the context of 
interference modeling from first principles for our domain of modelling, based purely on simple 
statistical assumptions on the system. 

B. Background on Univariate and multivariate a-Stable distributions 

In this section, we provide a brief technical survey of relevant results from the probability 
and statistics literature relating to sub-exponential family models of distributions, in particular 
the family of a-Stable models. These results will provide sufficient coverage to understand the 
derivations and results we develop in this paper related to interference modeling in wireless 
communications. 

Considered as generalizations of the Gaussian distribution, a-Stable models are defined as 
the class of location-scale distributions which are closed under convolutions. In an interference 
modeling context, a-stable distributions possess several useful properties, including the possibil- 
ity of incorporating infinite mean and infinite variance, skewness and heavy tails, see [21] and 
[20]. It is due to this inherent flexibility that they have found extensive use in practical modeling 
settings, both in wireless communications such as interference modelling that we consider and 
in many other domains of application, see a comprehensive list of such literature in the stable 
bibliography^. 

1) Univariate a-Stable Models: We consider a random variable X with a-stable distribution, 
denoted by X ~ S a (x; f3, 7, 5, 0). Where, S a (x; /3, 7, 5, 0) denotes the univariate four parameter 
stable distribution family under parameterization S(0) as defined in [22]. 

The univariate a-stable distribution we consider is specified by four parameters: a e (0, 2] 
determining the rate of tail decay; (3 e [—1,1] determining the degree and sign of asymmetry 
(skewness); 7 > the scale (under some parameterizations); and 5 E R the location. The parameter 
a is termed the characteristic exponent, with small and large a implying heavy and light tails 
respectively. In general a-stable models admit no closed-form expression for the density which 
can be analytically evaluated point-wise, except Gaussian (a — 2,/3 — 0), Cauchy (a = = 0) 
and Levy (a = 0.5, — 1) distribution cases. Therefore, statistical inference typically proceeds via 
the characteristic function, see discussions in [23], [24] and [25]. However, intractable to evaluate 
point-wise, importantly for wireless communication applications, simulation of random variates 
is very efficient (see [26]). 

Definition 1: A random variable X is stable if and only if X = aZ+b, where < a < 2, — 1 < /3 < 1, 
a > 0, b G R and Z is a random variable with characteristic function 




(1) 



1 http : // academic2 . american .edu/ jpnolan/ stable/ stable. html 
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where sign{u) = — 1 if u < 0, sign{u) = if u = and sign{u) = 1 if u > 0. 



From this definition one may define several practically useful reparameterizations, in this 
paper we consider the following parameterization denoted in [22] as the S(0) parameterization. 
A random variable X is said to have a stable distribution, S a (j3, 7, 5; 0), if its CF has the following 
form: 



E[exp(i9X)\ 



exp{-7 Q |#| Q (l + z/3(sign(0)) tan(^)(|7^| 1 - Q - 1)) + iSO) if a ^ 1 
exp{- 7 |0|(l + i/3(f)(sign(0))ln(7|0|)) + i^} if « = 1. 



In the following lemmas we present some fundamental basic facts about univariate a-Stable 
random variables that will be required to establish the results we develop in this paper. In 
particular these results will be used to construct analytic exact Poisson and doubly stochastic 
Poisson mixture representations of wireless communications interference processes, arising in 
the context in which an unknown number of interferers are present. This will be achieved by 
considering an important sub-family of a-stable models, those that are symmetric and isotropic. 

Lemma 1: If Y ~ S a (/3 : 7, 5; 0), then for any o ^ 0,6 6 R, the transformation Z = aY + b is a 
scaled version of the a-stable distribution. That is Z ~ S a (sign(a)(3, |a|7, ad + b; 0). In addition, the CF, 
densities and distribution functions are jointly continuous in all four parameters (a,/?, 7, 8) and in x. 
These results follow from [20] and [22, Proposition 1.16]. 

Lemma 2: If for all i G {1, ...,N} one has random variables Xi ~ S a (/3i, ji, 8f, 0) then the 
distribution of the linear combination, given N, is 

N 

Z = J2^~S(aJ,^,5-0) 



i=l 



N yiV n a 

7" = !>.". ?=4f$% (2) 

i=l L^ii=\ li 

Ef =1 8> + tan f (py - Ef =1 Pfli) { f « ± 1 > 

Eii ^ + f (/37 log 7 - E 1 1 lo § 7i) = 
This result follows from [20, Section 1.2, Property 1.2.1] and [22, Proposition 1.17]. 

A practically relevant sub-family of a-stable distributions is obtained when one considers 
the symmetric case. A random variable X is said to be distributed from a symmetric ct-Stable 
distribution, X ~ S a (0,7,5), when the skewness parameter j3 = 0. In this case, the model still 
captures a spectrum of distributions ranging from Gaussian a = 2 through to infinite mean and 
infinite variance models. This particular sub-class is interesting as it can be represented uniquely 
by a Scaled Mixture of Normals (SMiN) representation as shown in Lemma |3] 
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Lemma 3: In [27, Equations (3) and (4), pp.2] it is shown that when X ~ <S a (0,7, 5) it may be 
represented exactly by the following SMiN representation through the introduction of an auxiliary random 
variable A, 

X\X ~ N(5,-yX) , (3) 

with auxiliary scale variable X ~ S a / 2 (0, 1, 1). 

This result will be used in our models as the basis for working with a large family of symmetric 
stable models that we derive. Clearly, this is advantageous as conditionally on A, one has a 
Gaussian distributed random variable X\X. However, general closed-form expansions of the 
probability distribution functions in terms of well-understood functions do not exist. However, 
since the work of [28] who showed that all continuous stable distributions can be written in terms 
of infinite series expansions of elementary functions, there has been specific examples created. 
There is the series expansion representations developed in [29] they transform the parameters 
and provide the stable densities in terms of Foxs H functions, see [30]. Other reparameterizations 
to obtain representations include [31] who obtain representations with respect to incomplete 
hypergeometric functions. So in general for symmetric stable settings one may choose between 
a SMiN representation or a series expansion. A summary of the series and integral expansions 
for a-stable models is given by H} see [21] for details. 

Lemma 4 (a-Stable Density and Distribution Representations): w.l.o.g. the density function of an 
a-Stable distribution (standardized such that 7 = 1 and 5 = 0) can be evaluated pointwise according to 



f x {x;a,P,lAS{0)) 



±Re { J °°exp (-itx - t a exp (-i^/3K(a))) dt) , if a ^ 1, 
^Re { f °° exp (—itx — 1 t — i/3t logt) dt} , if a = 1. 

Alternatively via the series expansions [21] [Equation 2.4.6, p. 89] 



(4) 



f x (x;a,/3,l,0;S(0)) 



OO , 



-lHi+i) 
T(n+1) 



sin (77,7173)2; 



M-l 



n ~ x nb n x n - x , 



if a > l,/3 e [-1, l],x e 
if a = 1,(3 e (0,1], x e M 
if a < l,/3 e [-1, l],x e 



(5) 



where the coefficients b n are given by 

b n 



T(n 



exp (— /3u\nu) u n sin (1 + /3)u 



7T 

2J 



du. 



(6) 



In addition, the distribution function of an a-Stable model can be evaluated pointwise according to 

F x (x;a,P,l,0;S(0))-- 



(7) 



C(a,9) + J^exp (^— X"- 1 U a (cp, 9) j dip, if a ^ 1 and x > 0, 
\ f\ x exp ^— exp ^— I J f7i(«y3, /?)J d<^, if a = 1, and > 0, 

otherwise in all other cases it suffices to utilise the duality principle of infinitely divisible stable distributions 
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which has the consequence that 



F x {-x- a,P,l, 0; S{0)) + F x {x; a, -(3, 1, 0; S{0)) = 1. (8) 
Note, the notation of [21] [page 74] is adopted above in which 



e(a) = sgn(l — a), K(a) = a — 1 + sgn(l — a 
K(a) ^ m _ 1 1 
a 

sin f a(v + #) \ 1_a cos | ((a — l)ip + a 



= 0-^, C(a,9) = 1--{1 + 9)(l + e{a)), (9) 



\ COS 7}(f J COS ^</3 

To proceed we demonstrate when each approach will be of utility in the setting of interference 
modelling in wireless communications. 

Definition 2: Distributions of the same type will be defined to be those that differ only in location 
and scale. 

Without loss of generality, consider X a ~ S a (0, 1, 0), symmetric scaled a-stable random vari- 
ables. Variables X will have distributions which are of the same type as random variable X 
given in Definition [fl In the following specific cases we may replace the infinite scaled mixture 
SMiN representation with analytic functions representations as alternative models to evaluate the 
density of a symmetric a-stable random variate: 

Lemma 5: In the symmetric stable case, under two particular choices of a 6 ||, §}, the distribution 
of the stable model can be analytically represented via known functions according to: 

1) Finite Mean Interference Model (a = 3/2): the Holtsmark density [32] is given for random 
variable denoted X 3 / 2 which has an analytic density which is represented according to hypergeometric 
distributions as follows, 

1 „ (h 11 1 1 5 2 2 x 6 \ x 2 „ /3 . 5 2 5 7 4 2V 

7T 



fx {x) = -r(5/3) 2 F 3 — , — ; -,-,-; -— - - — 3 F 4 -, 1, A 



12' 12' 3' 2' 6' 3 6 J 3tt "V 4 ' '4'3'6'6'3' 

7z* . , „ /13 19 7 3 5 2 2 x 6 \ 
+ ^ r (4/3) 2 F 3 (-, -; -, -, -; -_ , -oo < x < oo, 



(10) 



with a = |,/3 = 0,7 = l,5 = and: 



p F q (ai, . . . , a p ; 6i, . . . , b q ; z) = ^ — 77^-7, (H) 

n=0 \ u l)n ■ ■ ■ \Uq)n It- 

with (a)„ = a(a + 1) . . . (a + n — 1) and (a) = 1. 
2) Infinite Mean Interference Model (a = 2/3): the Whittaker function density representation [33] 
is given for a random variable denoted X 2 / 3 which has an anlaytic denisty which is represented 
according to Whittaker functions, see [34] as follows, 

^ X ^ = 2vW' X ' _1 6XP i^7 x 2 ) W ~ 1/2 ' 1/6 {^f X J ' ~°° < X < 00 ' ^ 



with a = |,/3 = 0, 7 = 1,<5 = and 

We can a/so express W\-^ (z) in terms of confluent hyper geometric functions: 

W X]I1 (z) = exp(-z/2)z^ +1 ^U(fi - A + 1/2, 1 + 2/x; z) (14) 

^(a. ^ A = J (1 ~ h b \^ M ^ b > z ) + ^?^-^ h M{a -b+l,2-b;z) (15) 
1 (a — o + 1) 1 (a) 

and 

M(a,b]z) = xF^a^z) (16) 

The combination of Lemma [Q Lemma |2] and Lemma |3] will be directly utilized in results 
obtained in Theorem [2] In addition the results in Lemma [5] are of relevance in the practical 
examples developed to illustrate the theoretical results obtained. 

2) Multivariate a-Stable Models: It is also relevant to consider some background on multivari- 
ate ct-stable random variables, in particular the symmetric isotropic bi-variate a-Stable random 
variable. Multivariate a-Stable models are covered in detail in [22], [35], [36], [37]. Here we 
first present the joint CF for an elliptically contoured multivariate a-Stable distribution, see [38, 
Proposition 2.5.8], 

(fx = E [ju> T X] = exp (- (w T Sw) a/2 + juj T S\ (17) 

for some positive definite matrix £ and translation vector 6. Furthermore, there is a generalization 
of the result of Lemma [3] for the multivariate settings given in Lemma [6] below. 

Lemma 6: According to [20, Proposition 2.5.8], consider G ~ N(0, E) such that Y e IR d and 
A ~ S a / 2 (l, 7, 0; 0) as an independent univariate a-Stable random variable with < a < 2. Then the 
transformed random vector X = yXY is a-Stable and elliptically contoured with CF, 

fx — E [j^X] = exp ( (2) " /2 sec (^) (u^Eo.) . (18) 

This substable family reduces in the isotropic (radially symmetric) case, in which S z's diagonal, to: 

= E [ju; r X] = exp (( 7o ) Q |u,| Q + ju, r <5) , (19) 

for scale parameter 7o > 0. Furthermore, in this case the spectral measure is a uniform distribution on the 
unit-sphere and we obtain, 

X ~ N(6, AS) . (20) 
In addition, to the result in Lemma [6] for the generalization of the SMiN reprsentation to 
isotropic symmetric multivariate a-Stable random vectors, it is also relevant to present the general 
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properties for the density and distribution functions. 

III. System Model and Assumptions 

In this section we complete the system model introduced in the introduction and detail the 
required statistical assumptions that will be used throughout the paper: 

1) We assume a wireless network with N(t) transmitters at time t, distributed on a region 
1Z with area A n , at locations indexed by £(£) = {LW(t)}. =1 . Furthermore, we assume 
that the number of transmitters varies over time and space stochastically according to an 
inhomogeneous spatial-temporal Poisson process with intensity parameter X(x,y,t). 

2) The i-th potential interferer transmits an i.i.d wide band signal, represented by: 

K 



where fk is the subcarrier frequency and Xj? the source symbol from interferer i on sub- 
carrier k. 

3) The bandwidth, as quantified by the number of frequency carriers K, will be considered as 
stochastic with a truncated Poisson distribution given by K ~ Pois (A^) VJf £ {1, . . . , K max }. 
For example, this is practical if one considers that for all potential wideband interferers we 
assume one of the following two scenarios: all interferers transmit in the same bandwidth, 
but this bandwidth is unknown to the receiver and modelled according to a truncated 
Poisson distribution given by K ~ Pois {\k) Vif £ {1, . . . , K max }; alternatively, all interferers 
utilize the same total bandwidth, but the frequencies occupied by any given user may not 
overlap, however, the total bandwidth per user is unknown to the receiver and modelled 
according to a truncated Poisson distribution given by K ~ Pois (Xk) VK £ {1, . . . , K max }. 
In deriving the theoretical results, the expressions we obtain are interpretable in either of 
these practical scenarios. 

4) The random distance of the i-th potential interferer from the receiver is denoted by Ri and 
is given by: 



where is a random location of the i-th potential interferer and l R is a known location of 
the receiver in region 1Z. In a first part of the paper we assume that the location of the i-th 
potential interferer for the case of Models 1 and 2 will be uniformly distributed in space, 
such that the circular interference domain is given by: 




(21) 



fc=i 



(22) 




(23) 



10 



where tt is the maximal distance in which an interfere can have a non-negligible contribu- 
tion to the interference. 

5) For the z-th potential interferer, the low pass representation of the channel experienced by 

(i) A^e? k 

the symbol Xi ' is given by k _ a/2 ■ The path loss experienced by the i-th potential interferer 

is given by R { 2 , where a is the attenuation coefficient, a deterministic and known parameter 
reflecting the physical environment in which transmission is occurring. A^e^* is a complex 
coefficient that contains the shadowing and multipath fading (the amplitude distribution is 
not important and the phase is uniformly distributed over [0, 2tt]). 

6) After the adapted filter at the receiver side, the resulting total interference is given by: 

N K K N K N 

Y = E E Afxfcf^f = y. E + i E E C < 24 > 

i=l fe=l fe=l i=l fe=l t=l 

where cj^ is a random variable resulting from the filtering adapted for subcarrier k and 
depends on the system parameters. 

IV. Analytic Distributional Results for the Total Interference in homogeneous 

PNSC(a) Interference Models 

In this section we present the distributional results for the total interference given in (|24)> . We 
first detail explicitly the results for Model I in which we consider a homogeneous spatial and 
temporal intensity pattern for the distribution of interferers in the plane. Furthermore, we also 
consider the second order homogeneity of the process in which the intensity of the utilization of 
frequencies and therefore the number of occupied frequencies by each of the potential interferers 
is temporally and spatially homogeneous. This is significantly extended in section IV-BI where 
we present abridged generalizations of these results for the inhomogeneous in time and space 
settings. 

The results we develop in this section include Theorem [TJ which proves that at a given frequency 
of transmission, the total interference at the receiver experienced by interference from an unknown 
number of randomly distributed interferers can be shown to be an isotropic a-stable model in C 2 . 
This result is known, however it is informative to the models and extensions to present clearly 
a complete derivation in Appendix QIIJ We then significantly extend this result to provide three 
possible equivalent analytic representations of the distribution and density functions of the total 
interference at the receiver. In Theorem |2] we detail the PNSC(ct) model for a wideband system in 
which an unknown number of interferers, distributed randomly in the plane and transmitting on 
an unknown number of carrier frequencies is considered. This is achieved through development 
of representations of the resulting process as a Cox process, which in such models produces a 
Poisson-SMiN mixture representation. We provide two special examples of members from the 
family which characterize two important cases of infinite and finite mean interference models. 
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A. Model I: distributional results for total interference under homogeneous Poisson field of interferers in 
a circular domain 

In this section we consider the first model defined according to the following assumptions on 
the intensity of interferers in the field. 

Model 1: In this case we assume the intensity parameter X(x,y,t) = A, in other words, the 
mean number of transmitters does not change over time or space with distribution 

P (N(t) =n)= (Ay V exp- AA ^ . (25) 
n\ 

Given the system model for the total interference we now present the main result, which 
is to derive novel representations of the density and distribution in closed form for the total 
interference in Equation ([24]). 

Given a homogeneous spatial Poisson process with intensity parameter A, to define the distri- 
bution for the stochastic number of interferers in a given region of space, 

A {A n ) = {x G M 2 : < r < r T ) . (26) 

In this setting we show the total interference has closed form analytic density represented ac- 
cording to the results derived in Theorem [1] 

We begin by presenting a distributional result for the fc-th transmission frequency of the 2-th 
potential user, based on previously derived results for narrow-band systems in [8]. In this case 
we can write the characteristic function based on in-phase and quadrature-phase components 
presented in (|2i~]) according to definition |3l 

Definition 3: The characteristic function of the interference at the k-th transmission frequency, 
from i-th potential interferer, is given by: 

p y (M yr = E y(M [exp (jJMy™ + ju™Y™)] . (27) 

I Q 



Using (f27)) we can define for the A>th transmission frequency the characteristic function for the 
total interference, given an unknown number of independent potential interferes N in Lemma |7l 

Lemma 7: The characteristic function of the total interference for the k-th transmission frequency, 
for a random number of N {An) potential interferers, is given by: 



(k) (k) 
f Y (k) Y {k) \ U I ' U Q 



E 



R,c fe ,A fe ,# fe ,7V 



N 



X 



5>< 

i=l 



- a/2 Afcf C os | $ 



I 3 

(0 



OJ 



arctan 



(fc)V 

Q J 



(AO 



Given the expression for the CF in Lemma [7] we marginalize over the unknown number of 
interferers in region A^ to obtain the result in Lemma [8j 
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Lemma 8: The characteristic function, marginalized with respect to the random number of unknown 
■potential interferers, N, at the k-th transmission frequency, is given by: 



(k) (k) 



X 



N 



(fc) 



N=0 



8=1 



exp I j a/ I wf } 



(fc) 



(28) 



Equation (|23|) gives the spatial distribution of the interferers, conditional on the number of 
interferers present at the k-th frequency. Furthermore, in many situations (asynchronism, no 
power control), the signal strengths (included in the term c k ) of the interferers can be assumed 
independent and identically distributed. Hence, the CF for the total interference at frequency k 
in Lemma [8] is expressed according to Lemma HJ 

Lemma 9: Under the assumption that, given N potential interferers in region An, the spatial 
distribution of the locations for each of the interferers is given by density in (|23l > we can express the 
CF for the total interference at the k-th frequency according to: 

l(>yW tY W (wf^Wg^J = log ((PyW-yW \ U^ k \u^ 



exp i^R-^ 2 A k c k ^{uf ) 



(*)V 



x cos $/, — arctan 



(fc)' 
Q 



UJ 



Proof See Appendix H ■ 

In Lemma [10] we re-express the argument of the expectation as a complex series expansion in 
terms of Bessel functions and then marginalize over the random variable $fc in the expectation 
operator. 

Lemma 10: The log CF representation of the total interference at the k-th frequency, after marginal- 
izing the random variable for $ fe , is given by: 



I (*0 (k) 



(29) 



Proof: See Appendix [Q] 



Without loss of generality, we assume that the receiver, from which we measure the distances 
of each potential interferer, lies at the center of region An. Utilizing Lemma [7] through to Lemma 
HQ] we can now state the following result in Theorem [1] for the CF of the total interference at the 
k-th frequency. 

Theorem 1: The log CF, for the total interference at the k-th frequency, for a random number of 
potential interferers N in a region An, expressed in Lemma [TOl can be re-expressed in the form of a CF 
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representing the family of isotropic bivarite a-stable distributions S(a, 0, 7, 5; 0), 

</y j(fc , y(fc) ( w « , W «) = -7 A/ ( < ; ) + ( < ) • (30) 



Proq/: See Appendix [TTT] ■ 

Remark 1: Theorem\l]presents the bivariate log characteristic function of the complex random variable 
y(k) ^^(Yy*^ + Yq ) in (l24l >. T/ze resulting characteristic function is a member of the elliptic 
family of stable distributions for all < a < 2 and 7 > 0. In this case we get a = - and 7 = 



AttE 



A consequence of this result is that we may represent the density and distribution functions for 
Y ^ as a real random vector with the first component corresponding to its real part and the second 
component its imaginary part. We can then exploit one of the three common representations of 
the density for an isotropic bivariate stable vector, given in Lemma ITlT 

Lemma 11: Consider a multivariate isotropic a-stable random vector Y of dimension d with scale 7 
and location 6. Its density can be represented as: 
Representation I - SMiN The density f Y (y) is given by: 

/yia (y) = N (y; 8, AS) , (31) 
with A ~ S a /2 (A; 1,7, 0; 0) which is the S parameterisation of Nolan [35]. 

Representation II - Projection: [35] For every vector u e M. d , the one-dimensional projection (u, Y) is 
a univariate a-stable symmetric RV with stability index a. As detailed in [22], [35], the projection onto 
vector u in the isotropic case is given by the stable univariate random vector: 

(u,Y)~S a (0, 7 (u),<f(u);0). (32) 

By Cramer-Wold these univariate projections characterize the joint distribution, where 7 (•) and 5 (•) are 
called projection parameter functions, see definitions in [20], [21] and [35, Section 2.1]. In the special case 
of the isotropic multi-variate a-stable model one gets Vu e R d the simplification 7 (u) = 7. 



Given (f3~l) and (|32|) , we derive the expression for the density and distribution function of the total 



interference for an unknown number of interferers with unknown bandwidth for a few possible 
scenarios relating to how each potential interferer utilises the carrier frequencies available. To 
achieve this we utilise one of the three possible representations discussed in Lemma [IT] This is 
complicated by the fact that we will be working with the Poisson model under a temporal Cox 
process setting for the bandwidth utilised by each of the potential interferers. We first need to 
derive expressions for the distribution and density functions for the fc-fold convolution integrals 
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obtained from linear combinations of random vector Y K = J2 k=1 Y^: 

fc* (y) = ff* (y) = (/n * /y» * ■ • ■ * M (y) = {f^* * / Y ) (y) 

= J f¥~ lh (y - r) jfy (r) dr, 
and the CDF is expressed according to: 

F^k (y) = if * (y) = (F Y1 * F Y2 * ■ • • * F YK ) (y) = (i^'^* * F Y ) (y) 
= | 4^* (y " r) /y (t) dr. 



(33) 



(34) 



Remark 2: Analytic solutions for these convolution integrals can be obtained since we are working 
with a special class ofRVs characterized by the stable laws. This is a key advantage of stable RVs detailed in 
the univariate case in Lemma\2\ Generalizing this result to the bivariate setting developed in Theorem^and 
representations in Lemma [II] is non-trivial. In order to do it, we work under the projection representation. 

In Theorem [2] we evaluate the distribution of the total interference in (|24"]) . 

Theorem 2: Compound Poisson PNSC(a) Interference Model 
The PNSC(a) Interference Model can be expressed by the distribution of the total interference (indexed 
by a), in the case in which the unknown bandwidth of each potential interferer K ~ Pois (Xk) , e 

{1, • • • , Kmax}/ 

K N K N 

fc=l n=l fc=l n=l 

can be obtained uniquely by projection onto the real axes (radial symmetry argument). Therefore, Vu e 1Z 2 
the density and distribution of the projection RV Y (u) = (Y, u), are represented by 



/y (u) (y) = c K ] mx Y,F(K = k) f^ K(u) (y) I (k < K max ) 



k=l 
K„ 



(36) 



k=l 

where in the model defined in Theorem [0 and using Lemma [2j the following holds: 



f3 k (u) = 8 k (u) = 0, Vu G R d and k e {1, . . . , K max } , 

r°° j, (t\ 

dx, 



lk (u) = = k l a and 1 = ^A k , Ck \(A k c k )^} \ 

i=i L 



4 

a> 



(37) 



fc=i 

ITzz's representation admits the following closed form expressions, for the distribution and density function, 
according to Lemma [3] groen by: 
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1) The conditional density of Y (u) |A e M z's gz'yen by the following finite mixture of Gaussians, 

Kmax yk ^ / 1 \ 

/n«>i« (v) = %L E ex P (- A ^) -# 27rv ^(^: exp 1 {k ^ Km ^ ' (38) 

where k ~ S a / 2 (0, 1, 1; 0) . 

2) TTze conditional distribution function of Y (u) \k e R z's gz'uen by the following finite mixture of 



Gaussians, 



max 



Fy (u)lK (y) = P (y (u) < y|«) = £ exp (-A*) ^$ ) 1 ( fc ^ *««) ■ < 39 ) 

' k =i kl \ v ik (u) k y 

Proo/ - the proof of (|38|) and (|39|) is obtained in two steps. First, we take the spherically 
symmetric bivariate stable distribution derived in Theorem |2] and utilize the closure under con- 
volution result presented in Lemma |2] for each Poisson mixture component to obtain the result in 
(|36"1> . Then we use the projection and SMiN representations in Lemma [IT] to obtain a parametric 
closed form representation which can be evaluated. ■ 

Example 1 (Distribution of the PNSC(a) interference model): We illustrate in Figure\l\the model 
with the probability density function when the mean interfering bandwidth increases. We select 
a useful bandwidth given by K max = 64, and we choose different mean occupied bandwidths 
for interferers and plot the corresponding probability density functions (Subplots a, b and c). The 
dispersion parameter is set to 7 = 1 for this illustration. 

When the bandwidth occupied by interferers increases, the importance of strong users increases as 
can be seen on figureUi subplots d and e. As a consequence the peak at the center (y = 0) is reduced. 
The impulsiveness of the interference is increased if interfering users (for instance secondary users 
in cognitive radio) are allowed to use more bandwidth. 

Remark 3: The result developed in Theorem [2] z's significant as it provides a closed form analytic 
expression for both the density and distribution functions of the total interference for an unknown number 
of users and unknown bandwidth which were derived as a compound process. It is well known that in 
general closed form expressions for the distribution or density functions are generally non-analytic for 
general compound processes. Typically, evaluation of such distributions can be calculated pointwise via 
computationally expensive approximations, in the univariate case, via either Monte Carlo approximation 
or Panjer recursions, see discussions in [9]. 

Next we present two examples for the total interference based on lemma [5] which provide 
results when the interference for each frequency in the bandwidth has either an infinite mean 
(case a = 2/3) or a finite mean (a = 3/2). We assume that the number of interferers per frequency 
carrier is Poisson distributed with an intensity A and that the total number of frequencies occupied 
by all users, denoted k E {1,2,..., K max } is truncated Poisson with intensity X K . 

Example 2 (Compound Poisson PNSC(3/2) Holtsmark Interference Model): For any projection 
u e M. d the resulting doubly stochastic Poisson-Stable process is comprised of each mixture compo- 
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Fig. 1. PNSC(a) model for a — 1.5 and a = 0.5 with a useful bandwith determined by the number of carrier frequencies 
K max = 64. The subplots a to c give the distribution of the number of subcarriers used by the interferers. The larger Afc, the more 
bandwidth is used. Then we represent the tails and the center of the interference distribution for a = 1.5 (subplots d and/) and 
for q = 0.5 (subplots e and g). 



nent as defined in (33 in which we consider the case: 

(40) 



Y(u)\K = k~ f Y(u)lK=k (y) = S 3/2 (y- 0, Vkj, 0; 



with 7 given in ((37[) . Now defining the transformed random variable for Y(u) given K = k accord- 
ing to Zk = ;^=- Y(u)\K = k and utilizing Lemma\]\and Lemma^one obtains the analytic result 
for the density of Z expressed as a Poisson weighted mixture of Holtsmark densities: 

h (z) = c Kmsx Y: exp (-A*) ^ |-r(5/3) 2 F 3 -; -, -, -; -_ j 

fc=i k v 7 (41) 

- z - f (I 1 5 I + p H.Z ? h -- 2 ^f 

3tt 3 4 V 4 4' 3' 6' 6' 3' 3 6 j 3% 1 7 j2 3 ^12' 12' 6' 2' 3' 3 6 

Example 3 (Compound Poisson PNSC(2/3) Whittaker Interference Model): 
For any projection u G IR d the resulting doubly stochastic Poisson-Stable process is comprised of 
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each mixture component as defined in ((36)) in which we consider the case: 

Y(u)\K = k~ f YMK=k (y) = S 2/3 (y- 0, v^ 7 , 0; o) (42) 



with 7 given in (|37|) . Now defining the transformed random variable for Y(u) given K = k ac- 



cording to Z k 



Y(u)\K = k and utilizing Lemma\l\and Lemma\5\one obtains the analytic 



\fk~f 

result for the density of Z expressed as a Poisson weighted mixture of densities constructed from 
Whittaker functions 

Kmax \ k 1 / 2 \ / 4 \ 

1* (*) = c a-L E ex P (- A ^) ex P \tj z ) w -^ 6 \jj z ) ' (43) 



with Wx-fj, (z) given in (|T3|> . 



Next we present a result for the univariate a-stable distribution, for an analytic representation 
of the tail probability of the stable distribution in Lemma H2l 

Lemma 12: Given a random variable X ~ S(a, /3, 7, 5; 0) then as x — y 00 one can write f/ze limiting 
tail distribution 

F(X > x) ~ 7°c Q (l + flsx^oo 

(44) 

/x(x|o!, /3,7, 5; 0) ~ a 7 a c a (l + as x -»• 00. 

where c Q = sin (^) This result follows from ([22],Theorem 1.12) 

We can utilise this result to obtain analytic expressions for the tail probability of the Poisson 
Process for the total interference defined in (|24|). This analytic representation, presented in The- 
orem |3l is particularly useful for evaluation of tail probabilities and tail expectations, as may be 
relevant for calculations of BER, capacity and outage probabilities, ido to fix 

Theorem 3: The representation of the tail probability of the total interference, in the case in which 
the unknown bandwidth of each potential interferer K ~ Pots (Xk) , V/<" 6 {1, ... , K max }, defined by the 
probability 

(K N K N \ 

k=l n=l k=l n=l / 

can be obtained uniquely by projection onto the real axes (radial symmetry argument). Therefore, Vu e M 2 , 
the tail of the density and tail probability of the projected RV Y (u) = (Y, u), are represented as Y (u) — > 00 
by, 

fY(u) (V) ~ c K \ mx £ exp (-Ax) tJ« (% (u)) Q c a [l + (p k („))] (y)- (Q+1) , (46) 
fc=i 

and we may analytically present the tail probability by the limiting survival function given as Y (u) j 00 
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by, 

K N K N 



p(y ( U ) > y) = p x; E ^ (fc,n) + j E E ,n) > 

\fc = l 71=1 fe = l 71=1 / 

CkL E ex P (- A ^) ^ (T* ( u ))° c ^ + ( u )) )»" 



(47) 



fe=i 

where m f /ze model defined in Theorem [0 and using Lemma |2l f/ze following holds: 

(3 k (u) = S k (u) = 0, Vu e R d , and k e {1, . . . , 7^ (u) = £* = i if = kT, 

7 = \^ Hk ,A k \{H k A k y] / °° 4 Y:tr P{K = k) and c a = sin (f ) ^. 

Proof: The proof of the result for the representation of the distribution and density of the tail 
of the total interference follows by taking the spherically symmetric bivariate stable distribution 
derived in Theorem [TJ and utilising the closure under convolution result presented in Lemma 
|2] for each Poisson mixture component. Next, the representation of the tail density in (j46]) and 
tail distribution (limiting survival function) in (|47|) are obtained using the projection and SMiN 
representations in Lemma (TT] to obtain a parametric closed form representation which, in the 
limit, allows us to apply Lemma | 



V. Analytic Distributional Results for the Total Interference in inhomogeneous 

PNSC(a) Interference Models 

In this section we consider a significant extension of the model presented in the previous 
section to the Doubly Stochastic Poisson-SMiN models. It can be considered as stochastic intensity 
models in the Cox process setting. The results in Theorem H] present analytic distribution and 
density functions for the total interference in the case in which one models an unknown number 
of potential interferers, distributed uniformly in the plane and each occupying an unknown 
bandwidth in which the mean bandwidth of the potential interferers is stochastic and given by 
a Poisson-Gamma conjugacy. We then extend to the case of Model II and Model III in Theorem 
|5l in which we consider that the intensity of the interferers in the plane and in time can be 
inhomogeneous. 

A. Model I: stochastic intensity for the bandwidth occupied by interferers. 

Theorem 4: [Doubly Stochastic Poisson-Gamma PNSC(a) Model] 

Consider the case in which the bandwidth occupied by each potential interferer, quantified by the number 
of carrier frequencies K, is a random variable from inhomogeneous Poisson process: K ~ Pois {\k) Vi£" G 
{1, . . . , Kmax}, where Xk ~ Qa(a,b). The distribution of the total interference follows a temporal Cox 
process (doubly stochastic Poisson process) represented by a compound process given in (|24]> . Under any 
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projection, (Vu e R 2 ), the density and distribution of the projection Y (u) = (Y, u), are represented by: 

Mu) (v) = <**L E { a-m VT+b) \T+b) Sa ( y; ~ Pk (u) ' % (u) ' ~ 5k (u) ; °) • (48) 

k=l 

where the parameters were defined in $37) except: 



TTzz's representation admits the following closed form expressions, for the distribution and density function, 
according to Lemma [H given by: 

1) The conditional density ofY (u) |A G R z's gz'mr by the following finite mixture of normals: 



, ^ (g + fc-i)! / 6 y/ 1 y^^__ / 1 



a 



l)!fc! Vl + 6/ \1 + V 27rv/7 fc (u)A V 2^ (u) A 



(50) 



w/zere A ~ S a / 2 (0, 1, 1; 0) . 
2) The conditional distribution function of Y (vl) |A G R z's gz'uen by the following finite mixture of 



normals, 



F Y{u)]x (y)=F(Y (u)<y|A) 



-^l^t^x^'K^)' <5i) 

iwtfi F y(u) | A (0) = P (y (u) = 1 A) = exp (-A). 

Proo/ - See Appendix HVl ■ 

Remark 4: The result developed in Theorem [2] z's significant as it provides a closed form analytic 
expression for both the density and distribution functions of the total interference for an unknown number 
of users and unknown bandwidth. It is well known that closed form expressions for the distribution or 
density functions of such doubly stochastic compound process are generally non-analytic. Typically, such 
distributions can be evaluated point wise via computationally expensive approximations, in the univariate 
case, via Monte Carlo approximations. 

As special cases of the general analytic results for the distribution of the total interference 
presented in Theorem HJ, we present in the following analytic channel model results for the infinite 
mean and the finite mean interference models in which the number of interferers in domain Ar 
is stochastic and Poisson distributed and furthermore the number of carrier frequencies they 
occupy is also stochastic and modeled according to a doubly stochastic Poisson-Gamma-Stable 
process (i.e. a Cox-Gamma-Stable Process - a special form of a renewal process): the number of 
interferers per frequency carrier is Poisson distributed with an intensity parameter A and that 
the total number of frequencies occupied by all users, denoted k G {1,2,..., K max } is truncated 
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Poisson with a stochastic intensity parameter \ K ~ Qa(a,b). 

Example 4 (Doubly Stochastic PNSC(2/3)-Holtsmark Interference): 
For any projection u £ R. d the resulting doubly stochastic Poisson-Stable process is comprised of 
each mixture component as defined in ([36]) in which we consider the case: 

Y(u)\K = k ~ f Y{u)]K=k (y) = S 2/3 (y- 0, v^ 7 , 0; o) (52) 

with 7 given in ^7) . Defining the transformed random variable for Y(u) given K = k according to 
Y(u)\K = k and utilizing Lemma\l\and Lemma\5\one obtains the analytic result for the 



density of Z expressed as a Poisson weighted mixture of Holtsmark densities 

iX (a + k- 1)1 ( b \ a ( 1 \ k f 1 ^ „ / 5 11 1 1 5 2V 



' (a-l)\k\ \l + bj \l + bj w '* \12' 12' 3' 2' 6' 3 6 

z 2 /3 5 2 5 7 4 2V\ 7 ^ m/ ox p f 13 19 7 3 5 2 ^ 6 
37 3 4 U' '4 ; 3'6'6'3 ; ~"36"J + 3% ( / )2 3 V 12' 12 ; 6' 2' 3 ; 3^" 



(53) 



Example 5 (Doubly Stochastic PNSC(3/2)-Whittaker Interference): 
For any projection u £ R. d the resulting doubly stochastic Poisson-Stable process is comprised of 
each mixture component as defined in ((36)) in which we consider the case: 



Y(u)\K = k ~ f Y (u)\K=k (y) = S3/2 (y; 0, Vk~-?, 0; 



0) (54) 



with 7 given in (f37)) . Defining the transformed random variable for Y(u) given K = k according 
to Zk = -7=- Y(u)\K = k and utilizing Lemma\l\and Lemma\5\one obtains the analytic result for 



Vk'y 

the density of Z expressed as a Poisson weighted mixture of densities constructed from Whittaker 
functions 



JzW-Ok^JL, ( a -l)\k\ [l + bj \l + b) 

leXp (^i^AVa(^- 2 ) 



(55) 



2V37T 



with Wx-fj, (z) given in (fl3|> . 



B. Models II & III: distributional results for total interference under inhomogeneous Poisson field of 
interferers in a circular domain 

In this section we generalize the model assumptions related to the functional form of the 
distribution of the interferers in time and space according to one of the following two model 
choices. 

Model 2: we assume the intensity parameter is X(x,y,t) = A(t). In other words, the mean 
number of transmitters is inhomogeneous in time and homogeneous in space with distribution: 

P ([N(t) - N(t - r)] =n)= V [*-*■■*] e -A [t _ T[f] Ag (56) 
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where X[ t - T ,t] = J t ^ T X(t)dt. 

Model 3: we assume the intensity parameter is X(x,y,t) = X(x,y), in other words, the mean 
number of transmitters is homogeneous in time and inhomogeneous in space. We consider the 
Borel o"-field of IR 2 in order to define a spatial Poisson inhomogeneous point process (SPIP). 
This allows one to define on a measurable space S C IR 2 , a random countable subset II. It 
will be governed by a stochastic mechanism that induces two properties for random variables 
N(A) = |n fi A\, i.e., the number of points of II lying in measurable subsets A of S. Specifically, 
for any finite collection A 1 ,...,A k of pairwise disjoint measurable subsets of S, the random 
variables N(Ai), . . . , N(A k ) are independent and for any measurable subset A of S, N(A) follows 
a Poisson distribution with mean J A A(x)dx. Here, A(x) is the intensity function of the spatial 
inhomogeneous Poisson process, a non-negative measurable function defined on S such that 
J A A(x)dx < 1 for all bounded subsets A of S. Thus the distributional properties of a SPIP are 
determined by its intensity function A(-), or, equivalently, by the mean measure of the process 
A(A), defined for all measurable subsets A of S. For theoretical background on spatial Poisson 
processes, see, for instance, [39], [40], [41]. This results in a SPIP distribution given by: 

P (N(A n ) =n)= ^ e- A ^. (57) 

n\ 

Given the system model for the total interference we now present the main result, which 
is to derive novel representations of the density and distribution in closed form for the total 
interference in (j24]) . 

Consider a setting in which we have a inhomogeneous temporal or spatial Poisson process 
with intensity parameter X(x,y,t) which defines the stochastic number of interferers in a given re- 
gion of space at a given time interval. In the first instance we consider the parametric interference 
space defined by: 

VI (A n ) = {x G M 2 : < r < r T } . (58) 

Under the assumptions of Model II and Model III, we generalize the analytic distributional 
results for the total interference developed for the homogeneous Poisson model derived in The- 
orem |4j We utilize a generic property of Poisson processes provided in Lemma [13] to re-obtain 
analogous distribution results under inhomogeneous Poisson models. 

Lemma 13: According to the mapping Theorem of [40, pp. 18], any inhomogeneous Poisson processes 
can be made homogeneous in space or time via a suitable monotonic transformation. Consider such 
inhomogeneous Poisson point process (II) in M. d x IR + with mean intensity measure A(-) defined over 
a spatial region A R eR d and time interval r G IR + with X (A R , r) = f T J A A(x, y, t) dt dx dy and define 
a mapping f :M. d x R + — y T. Given a smooth bijective function f then the transformation of the Poisson 
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process f (II) is itself a Poisson process on T with induced intensity measure X* given for any B C T 

A* (B) = A (r 1 (B)) . (59) 
Under Model II in which the intensity function is time inhomogeneous as specified by A (x, y, t) = 
X (t), one may utilise the result in Lemma [13] to transform the unknown number of interferers 
which is time varying into a time homogeneous Poisson process. The same is true for Model III 
in which the intensity function is space inhomogeneous as specified by A (x,y,t) = X (x,y). We 
illustrate this concept for Model III with the results presented in the following Theorem: 

Theorem 5: Consider Model III in which the unknown number of potential interferers is a spatially 
inhomogeneous Poisson process with spatial intensity function given by one of the following two scenarios: 
1) Scenario I: consider a mean spatial intensity for the number of potential interferers that decays 
according to a power law as a function of distance from the base-station. Therefore working in 
polar coordinates with X (r, 0) = Xqt 13 ^, for some coefficient of decay {3, base intensity measure A 
and r > . Defining u = r 13 we have the mapping from polar to Cartesian coordinates given by 



x = lo 1 ^ cos ((f)) and y = oo 1 ^ sin (0). Therefore when the inhomogeneous Poisson process, with 
intensity function X (r), is mapped to the plane (u, 4>), we obtain the intensity measure in the plane 
given by: 

X* = J Jx{r)dxdy = J J -r^oT^do; d<f>. (60) 

Therefore ignoring <p the homogeneous Poisson process r 13 has rate X* = A {^j 
2) Scenario II: consider situations where we have directional antennas. In this case we work with a 
restricted portion (sector) of the plane in which potential interferers may be present. Furthermore, 
assume that the interferers are Poisson distributed with mean intensity X restricted to this sector. 
Working in polar coordinates we consider this sector of the plane traced out by angle <fi, and we map 
the Poisson process in this sector to homogeneous Poisson process on the whole plane with rate: 

= (61) 



2ir J T J \ 27r . 

Therefore, working under this transformed Poisson process for the potential number of interfrers in the 
plane, we obtain the log CP for the total interference at the k-th frequency, given by: 



,/. (, .(*) , 



Proof: The proof of this result proceeds directly according to the result in Theorem |ITl] 
Corollary 1: According to Theorem |5l the results stated in Theorems |2] and |4] can be eas- 
ily derived under the transformed homogeneous Poisson process for the number of potential 
interferers in the plane. 



(62) 
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VI. Generalised SNR in PNSC(a) Interference Models 

In this section we consider an extension to the model proposed in [42] which develops the 
concept of an Additive White Symmetric a-stable Noise (AWSaSN) model. In this section we 
consider extending the analyse they performed to the generalized model developed in this paper. 
We now allow for an unknown number of potential interferers who are each transmitting on an 
unknown random bandwidth and interterfering from uniformly and unknown locations in the 
plane as specified in the results developed in Theorem [T] and Theorem |2l Again throughout this 
section we consider the results in the case of projection Y(u), which hold for all projections u e C. 

A. Generalised Geometric SNR in PNSC(a) Interference Models 

To begin the definition of the Geometric Signal-to-Noise Ratio (GSNR) for the PNSC(ct) interfer- 
ence model in a wideband system. We first define the following properties of a stable distributed 
random variable Y known as the fractional moments as defined for symmetric ct-stable models 
as follows. 

Definition 4: a-Stable Fractional Lower Order Moments (S(O)-FLOM's): 

Given a sub-exponential a-Stable distribution Y ~ S a (y; 0, 7, 5; S(0)) then the FLOM can be evaluated 
analytically according to the following moments or Fractional Lower Order Moment conditions for p e 

(-1,2), 



+00, p > a 



E[\Y\ 



2T(l+p)smipf™(l-Re{$(t)})tP~ 1 dt, O<p<a<l,0^O (63) 
C(p, a)j*, -1 < p < a < 2, = 0, 

where $(t) represents the characteristic function of r.v. Y and 

g fe.) = ^' r ffi r <- 5) , 

«V7rl (-|J 

and the following special cases of the first integer moment given by 

f<5-/37tan^, p = 1, 1 < a < 2 
E [Y] = I 2 (64) 

I 5, p — 1, a — 2. 

An alternative representation of the FLOM's results are derived for the parameterization S(l) 

in [20] and gives the results for the FLOM's according to Lemma [5] 

Definition 5: a-Stable Fractional Lower Order Moments (S(l)-FLOM's): 

Given a sub-exponential a-Stable distribution 

X ~ S a (x; 0, 7, 5; S(l)) in which a G (0, 2) and = when a — 1. Then the FLOM can be evaluated 
analytically according to the following moments or Fractional Lower Order Moment conditions for p e 
(0, a) according to 

E[\X\ p ]=C(p,a,0;S(l)f 1 (65) 
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where [20, Property 1.2.17] showed that 

C (p, a, (3; S(l)) = — ^ \ - ° J — ( 1 + (3 2 tan 2 —) 2a cos (- arctan ((3 tan —)) . (66) 

p J Q u _p_1 sin udu V 2 / Va V 2 // 

Therefore we can derive the following Theorem [6] for the FLOM's of the PNSC(a) Total 
interference model in three cases, the homogeneous case, the doubly stochastic Poisson-Gamma- 
PNSC(ct) interference model and the general inhomogeneous in time or space Inhomogeneous- 
Poisson-PNSC(a) model. 

Theorem 6: The FLOM's of the total interference in the PNSC(a) homogeneous interference model 
in Theorem [2] are given by 

{oo p > a, 

k / ^^ ( 67 ) 
C(P, a)c- K l ax Ef =1 exp (-X K ) ^ (Vk^ a < p < a. 

In the case of the doubly stochastic Poisson-Gamma-PNSC(a) interference model the FLOM's are given 
by 

{oo p > a, 

N £ (68) 

In the case of the inhomogenous in space or time PNSC(a) interference models, replace \k with X* in 
Theorem [13] 

Proof: These results follow trivially from application of the mixture representations of the 
interference models. ■ 

Having derived the FLOM's for the PNSC(a) Total interference models in wideband transmis- 
sion, we now consider the definition of the Geometric SNR in the PNSC(a) model. As discussed 
in [42] in an a-stable based interference model, the notion of noise power is mathematically 
undefined and the standard signal-to-noise ratio (SNR) is not a well-defined strength measure. 

The solution proposed involves a new indicator of the process strength, which is the geometric 
noise power. Since this estimator is a scale parameter it can be used as an indicator of process 
strength or power in situations where second-order methods are inadequate. The power is defined 
in [42, Equation 2] to be given in the total interference model for the PNSC(a) interference models 
according to Proposition [61 

Definition 6: For the interference Y the general definition of the geometric noise power which provides 
a measure of process strength is given by 

SO = S0(Y) = exp (E [log \Y\\) . (69) 
We note that using the FLOM's an upper bound on the noise power can be obtained via 
Jensen's inequality as given in Lemma H4l 

Lemma 14: Consider Y as the total interference in the homogeneous PNSC(a) or doubly stochastic 
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Poisson-Gamma-PNSC(a) models. Then the following bound on the noise power applies for p = 1 and 

a e (1,2), 

SO<E[\Y\]. (70) 

with E[|y|] as given in Theorem® 

In addition, we can also evaluate the noise power when Y is an a-stable random variable from 
the PNSC(ct) homogeneous total interference model or the doubly stochastic Poisson-Gamma- 
PNSC(a) model as shown in Lemma H5l 

Lemma 15: Consider Y as the total interference in the homogeneous PNSC(a) or doubly stochastic 
Poisson-Gamma-PNSC(a) models. Then the following evaluation on the geometric noise power applies 

so = %L Ef=i ex P (-**) ^Vk-fC* (7i) 

where C g is the exponential of the Euler constant. In the case that we consider Y as an a-stable RV from 
the doubly stochastic Poisson-Gamma-PNSC(a) homogeneous total interference model the total geometric 
power is given by 

g Q _ d K max Efc=l ( ( a -"i) !fc ! (iTb) (lTb) ^1 C 9 ^ 

In the case of the inhomogenous in space or time PNSC(a) interference models, replace X K with X* in 
Theorem [131 

We can then define the GNSR for each of the PNSC(a) total interference models in the 
wideband system according to the definition in [42] [Equation 4] is given in Proposition [161 
Lemma 16: For the interference Y the general definition is given by 

GNSR ^i0 (73) 

Therefore when Y is an a-stable RV from the PNSC(a) homogeneous total interference model the total 
GNSR is given by 

GNSR =^\-, k * r- A ■ (74) 

where A is the amplitude of a modulated signal. In the case that we consider Y as an a-stable RV from 
the doubly stochastic Poisson-Gamma-PNSC(a) homogeneous total interference model the total geometric 
power is given by 

GN SR = -L- | *Cg j (75) 

\ a K max 2^k=l (o-l)!fc! \l+b) \l+b) VK^g f 

In the case of the inhomogenous in space or time PNSC(a) interference models, replace Xk with X* in 
Theorem f 
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Remark 5: As noted in [42] the normalizing constant 2C 9 ensures that for the Gaussian case (a = 2), 
the GSNR becomes the standard SNRfor the PNSC(a) model 

In Example [6] we present a study of the characteristic of the GNSR for the total interference 
under the PNSC(a) model developed. This is important to analyze, since the properties of the 
a-stable model are known to be complex functions of a, (3, 7 and K. In this case the a-stable 
model is symmetric removing (3 = 0. 

Example 6: In the following simulations we consider the PNSC(a) homogeneous model with 
a range of a G [0.1, 1.9] values including three of which correspond to the special cases of the 
PNSC(a) Interference model developed which are the Compound Poisson PNSC(3/2)-Holtsmark 
Interference Model, the Compound Poisson PNSC(1)-Cauchy Interfence Model and the Com- 
pound Poisson PNSC(2/3)-Whittaker Interference Model. We set the GNSR according to the 
scale parameter given by 7 = {0.1,10,250,500,750,1000} x 10~ 5 , A = 1, K = 64 (64 carriers) 
and \k = 10. The results are provided in Figures |2] which demonstrate that there is a clearly non- 
linear relationship between a and the GSNR, with a linear relationship between 7 and GNSR as 
expected from the derive GNSR in Lemma |T6l 



GSNR for Wideband PNSC(cc) Interference Models vs. Heavyness of Impulsive Noise 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



a-Stable Tail Index 



Fig. 2. Relationship between GSNR,Tail Index a and scale 7 in the Compound Poisson PNSC Interference Model. 



VII. Likelihood Ratio in additive PNSC(a) interference 

We consider that the input X comes from a BPSK modulation and is affected by an impulsive 
noise Y modeled by the results developed in Theorems [T] and |2j Furthermore, we develop results 
below conditional on a fixed bandwidth, ie. fixed k, in the transmission and we assume that the 
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user transmits in the presence of an unknown number of potential interferers who are uniformly 
distributed in the field of transmission. 

In this scenario we consider two possible hypothesis for the received sample k, denoted by 
R k = X k + Y k with Y k = YliLi Yj k '^ + j YliLi Yq ,% \ given, according to Lemma [T] and Theorem S 
by the following: 

' Ho : R k = X(H ) + Y k ~S a (0,7 fc , 1; 5(0)) 

(76) 

Ux : R k = XiUx) + Y k ~S a (0, 7 fc , -1; S(0)) . 

Based on the previously presented representation of the distribution, we give analytical ex- 
pressions for the likelihood ratio in Theorem 

Theorem 7: Consider the BPSK transmitted signal in the presence of the homogeneous PNSC(a) 
interference model. The log likelihood ratio A(r) is given by: 

\( r \ - S a (0,7 k ,x(Ho);S(0)) 

which can be expressed according to the following forms depending on the tail index a: 

1) As the tail exponent of the stable distribution approaches 1, the distribution of the test statistic for 
any value of received signal reR converges in distribution to the Cauchy receiver with test statistic 
of [3 ], given by 

v am , Ik + {r - x{H )f 

km A(r) ->■ ^— -, tttt^- ( 7 §) 

et-H 7 fc + (r - x{Ui)Y 

2) In the infinite mean interference models characterized by the Compound Poisson Holtsmark Inter- 
ference Model, under an affine transformation of the received signal Z k = R k for the null and 
alternative hypothesis, the test statistic for the LET is analytically evaluated according to 



Um Mz) ^ ^(5/3)F 2 , 3 (z - x(Hq)) - ^g^F 3 , 4 (z - s(K )) 
o^a/s ir(5/3)F 2 , 3 (z - x(Hi)) - ( ^ 3 ( f ))2 F 3A (z - x(H x )) 

+r(4/3) 7( %y o))4 F 2 ,3 (z-x(Hq)) 
+r(4/3) 7(2 - 3 ^ l))4 F 2 , 3 (z - 



(79) 



with x(U.) = \-^\x(U.), x{U.) G {-1,1}, F 2t3 (z) = 2F3 |, if;-^) and F 3A (z) = 

P ( 1 1 5.2 5 7 4. 2 2 z 6 \ 
3 r 4 J-, 4 , 3 , 6 , g, 3 , 3 6 J- 

3) In the finite mean interference models characterized by the Compound Poisson PNSC(a) Interference 
Channel Model, the test statistic for the LRT is analytically evaluated according to: 

lim A(r) -> k ~ ^)|- 1 e^C^C^o))-^_ 1/a;1/e ( ^ (t . _ a? (^))~ a ) _ 
a ^ 3 / 2 \r - x{U{)\- 1 e& r - x ^ W„ 1/2;1/6 (±{r-x(Hi))- 2 ) 

4) In the case that one considers a general a e [0, 2] then the test statistic for the LRT after standardiza- 
tion involving scaling by ^~ 1 / a is evaluated (via truncation of sums) according to series expansions 
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provided in [21, Sections 2.4.6, 2.4.8], 



A(r) 



as— 1 



(r-x(n )) 2 + l 



-as—l ' 



< a < 1, 



(r 

E 



oo 
s=0 



x(Hi)Y + l 
D 1 / (r - x(H )) 2s 



a 



E^'(r-s(«i)) 

exp , _ 



2s ' 



1 < a < 2, 



(81) 



cxp 



a 



2. 



with D[ 



-T(as+l)sm(2f-), D 



ii 



(-!)« 



(2s)! 



Proof: The analytic LRT representations are derived using results from theorem |2] and theorem 
[T] for BPSK modulation. In addition, for item 4, the result follows from utilization of a series 
expansion representation of the a-stable model. ■ 

In the following simulation we consider the PNSC(a) homogeneous model. We evaluate the 
log likelihood ratio (LRT) at a range of values reR versus tail index a under two settings. The 
first is an exhaustive simulation via Monte Carlo methods, which indicates the true estimated 
LRT. Then the analytic expressions derived for different a is done. The limitations come for the 
Compound Poisson Holtsmark Interference Model f79~l> in the evaluation of functions p F q (.) in 
(ITT) and in the case of a general a (|8~TT> due to the sum truncation. 

Example 7 (Log likelihood ratio for a = 1 (fTS")) ): In this first example we plot the LRT in the 
well known case a = 1 (Figure [3]). As expected the fit is very good. This plot essentially allows to 



ct=l (Cauchy density) 




-6 -3 3 6 

r 

Fig. 3. Comparison between exact expression {78) of LRT versus Monte Carlo under the Cauchy model with a = 1. 

validate the accuracy of the Mote Carlo LRT evaluation. 

Example 8 (Log likelihood ratio for a = 3/2 (|79|) ): We plot the LRT in the case a = 3/2 (figure 
0) corresponding to the finite mean case - Holtsmark density. The fit is very accurate in the range 
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ct=1.5 (Holtsmark density) 




-6 -4 -2 2 4 6 



Fig. 4. Comparison between exact expressions of LRT versus Monte Carlo under the PNSC(3/2) model. 



[—3, 3] but for larger values the calculation of the functions p F q in ([79]) can not be performed accu- 
rately. It can however be sufficient in practice because larger values will be clipped at least due to 
the non linearity of the circuits. 

Example 9 (Log likelihood ratio for a = 2/3 (|80~1) ): We plot the LRT in the case a = 2/3 (Figure 
O corresponding to the infinite mean case - Whittaker density. The fit is very close except in very 



a=2/3 (Whittaker density) 




Fig. 5. Comparison between exact expressions of LRT versus Monte Carlo under the PNSC(2/3) model. 

small intervals around 1 and — 1. With a minimum care, this limitation is not a strong difficulty. 
Similarly to the Holtsmark density this distribution can give very interesting insights in the impact 
of impulsive noise on communications and allow to derive appropriate solutions for detection and 
estimation problems. 

Example 10 (Log likelihood ratio for 1 < a < 2 (181) ): We plot the LRT in the cases a = 1.8 and 
a = 1.2 (Figure® from the general series expansion case in (|81) . We see a very accurate fit between 
the real LRT and the truncated sum. The range of validity however depends on the number of terms 
we consider in the sum. For low values of r it converges very quickly (20 terms are sufficient) but 
we reach the computer limitations for larger values of r if we try to compute more than 400 terms 
in the sum. The range of convergence is also dependent on the value of a. The larger a is, the larger 
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Fig. 6. Comparison between an estimated LRT from truncated series expansions versus (Monte Carlo - exact) under the PNSC(a) 
model with a — 1.9 (left subplot) and a = 1.4 (right subplot). 



the range. When getting close to one, the series expansion (due to computer limitation) can only be 
used in a range smaller than [—1,1]. It can however have interest for instance for detection of small 
signals in a stable noise. 

Example 11 (Log likelihood ratio for < a < 1 (181) ): We finally plot the LRT in the cases a = 
0.8 and a = 0.2 (figure 0) from the general series expansion case in (|8T|) . Again we see a very good 

01=0.8 0=Q2 




-4 -2 2 4 



Fig. 7. Comparison between an estimated LRT from truncated series expansions versus (Monte Carlo - exact) under the PNSC(q) 
model with a = 1.2 (left subplot) and a = 0.2 (right subplot). 

fit between the real LRT and the truncated sum except for a small intervals around 1 and —1. The 
range where the model is not accurate gets larger as a gets closer to 1. However it remains quite 
small and, out of this range, less than 20 terms in the sum give an accurate result. 

VIII. Conclusion 

This paper deals with interference model that can be useful for ad hoc networks or cognitive 
radio. We consider a Poisson field of interferers with varying occupied bandwidth. This results 
in a doubly stochastic Poisson model. The resulting law is an a-stable distribution. Our main 
contributions are as follows: 
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1) the extension of the model with a random number of interferers to a doubly stochastic 
model; this includes a system with a random number of interferers having a random 
transmission bandwidth. It can have a direct application to cognitive radio systems where 
interferers are secondary users that adapt their bit rate to the environment they sense and 
where the primary users use an OFDM transmission like in LTE or IEEE 802.11; 

2) proposals for analytical expression of the probability density function and the cumulative 
distribution based on series expansion. Such tools can be very useful when tackling other 
problems like the capacity evaluation of large networks or solutions for detection and esti- 
mation. Our approach is based on two representations of the stable distributions: the SMiN 
representation, where the stable random variable is represented by a normal distribution 
with a random variance following a stable law, and the projection of a stable random vector 
which is a univariate stable random variable. 

3) we also include some special cases that can be obtained without the need to truncate the 
infinite sum of the series extension in two cases (other than Gauss and Cauchy which have 
exact expressions). We include a finite mean case (a = 3/2) based on Holtsmark distribution 
and an infinite mean case (a = 2/3) using Whittaker functions. Those expressions can be 
very useful in the sense they represent significant situations that could be used as reference 
for any evaluation of the considered scenario; 

4) the model is finally extended for inhomogeneous in time or space Poisson process. This 
allows to consider varying mean bandwidth occupation for interferers or spatially inhomo- 
geneous positions of interferers. 

We give an illustration of the proposed representation through the evaluation of the likelihood 
ratio for a binary source in additive symmetric a-stable interference, highlighting the accuracy 
of the approach, the rate of convergence of the involved series expansion but also its limits. 

Our results can have significant application linked to information theory, digital communica- 
tions or signal processing when the noise has an impulsive nature. 

Appendix I 
Characteristic function 

Proof: Under the assumption that, given iV potential interferers in region A n , the interferer 
locations are uniformly distributed, we can express the CF for the total interference at the /c-th 
frequency according to, 



N=0 



x cos $i. — arctan 




IV 



F(N(A n )) 
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Using the Taylor series representation of an exponent, we rewrite the CF in the log domain, 
according to: 



Avrr T E RiCfciAfei# , 



exp ^'i?- CT /M feCfc ^f 



(*) 



x cos §u — arctan 



UJ 



(ky 



UJ 



(fc) 



-1 . 



Appendix II 
Proof of lemma 8 

Proof: First re-express the argument of the expectation in the log CF in Lemma [8] using the 
complex series expansion based on Bessel functions, given by [34] 

oo 

exp (ja cos (9)) = '^^j s e s J s (a) cos (s9) 

8=1 

where eo = 1 and e s = 2 for all s > 1, and J s is the Bessel function of order s defined by: 

i r 

J s (x) — — / exp (— j (st — x sin (r))) dr. 
2vr J-n 



Applying this identity allows us to re-express the argument of the expectation in (|29|) according 
to, 



exp I jR zA k Ck\l [wf^ ) + i^n' ) cos I — arctan 



,(*) 



UJ 



(ky 
Q 



UJ 



(k) 



s=0 



^2j s e s J s I R *A k c k J + ( w Q fc) ) ) cos ( s ® k ~~ sarctan 



UJ 



(ky 
Q 



UJ 



(*) 



which allows us to write the log CF as: 

ijj (k) (fc) [UJ^.U^ 



x cos s<& k — s arctan 



ijj 



(ky 



UJ 



(k) 



- 1 



Next, under model assumptions in Section [TIT] we have that, for the A;-th frequency, the random 
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variable is uniformly distributed in [0,27r]. Therefore: 

/ oo 



,s=0 



j s ^/ 2 A fcCfe J(4 fe) ) +(4 fe) 



xE$ 



cos s$t — s arctan 



Co' 



(ft)' 

Q 



(ft) 



which holds due to the observation that E^ 



for s > 1. 



Appendix III 
Proof of Theorem I 

Proof: To derive this result for the CF as a representation of a isotropic bivariate a-stable 
distribution we make use of [43, Identity 3.12, p. 152] given by: 



lim a 



a 2r 

— exp (julOr)) dr — 1 



(/ 1 (x)) 2 juexp (jux) dx 



(82) 



Conditional on any given number of potential interferers, in the spatial region An, we can 
marginalize the CF given in Lemma [lOl with respect to the unknown spatial locations of these 
interferers. To achieve this we utilise the assumption on the spatial distribution of these interferers 
given in model assumptions in Section [Hi] Hence, we integrate the log CF given in Lemma 
as follows, 



, (ft) (k) 



Xnr'i 



I'T 



E, 



2r , 

-rrdr - 1 



T 



I'T 



^ E c fc ,A fc 


Jo ( 


\®c k ,A k 


J ( 



r 



-a/2 A 



k c k \j \ uf ] 



ijj 



(fc) 



T 



2r , 
— dr 



r 



T 
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Next we integrate by parts as follows, 



cir k ' k 



Avrr T E CfcjAfc 



J | r T a/2 A k c k \j ( uj 



dr Ck ' Ak 



Jo r-^ 2 A 



kCk 



J I r T a/2 A k c k 



- 1 



dr 



tt 2r 
—dr 

r T 



dr . 



To proceed, we next expand the region in which the interferers are distributed via the limit 
r T — > oo and apply the aforementioned identity 



lim \7ir T E Cfc! A, 



J I r T a/2 A k c k \ [uf ) 



o d r 



-E 



c fc) A fc 



J, (r-/»^( W ?>) , + ( W «) ! 



dr 



lim A7rr r E CfeiAfe 



J [ r T a/2 A k c k \l (w} K M + (w, 



lim A7rr 



r 



I'T 



dr 



E, 



J (r-/>^( W ?')'+(^') : 



dr. 



We can now evaluate the limits in each term above. Starting with the first one, we utilize the 
result from [1, Equation (12)] which allows us to state the following equivalent limit expression 
for the characteristic function for the total interference at two extremes: 



lim r T E c Afe 



Jn 



lim r T E CijAfe 

tt—>-0 



Now we note that since lim r% = and 

r T ->n 



limE Cfc)A)t 



r T ->-0 



Jol4 /2 A fcCfc W(<M +f w yj 1-1 



,(*) 



,( fc )Y 



0, we can apply L'Hopitals rule and the 
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identity A.j (x) = —J\ (x) given in [34] in conjunction with the chain rule to obtain 



lim r T 2 E C( , A( 

ry— »0 



——2 

= lim crr^ E c , a 

». 1 K ' 



—Jx [4A kCk J (u^y + (^ fe) 



A k c k \\ 



(*) 



(fc) 



We note that for a > 2 this limit converges to 0. Working with the second term, we have the 
following 



r r d 

lim A7r?4 / -r^c k ,Ak 



(*0 



= lim A7irf. 



E 



Jo r-^cJ^f)' 



2 



dr. 



Using the specified identity and noting the result of [44, Equation 17, pp. 8], we obtain 



, / (k) <k)\ 

V Yl (k) ,Y^ \ u i ,U Q J 



= -Att 



E 



4 



Ji (x) 



(83) 



\ v y v y / L J Jo 0> 

Equation (|83|) is the log characteristic function of an isotropic bivariate symmetric a-stable 
distribution, where the characteristic exponent a — -, and the dispersion parameter 



7 = A7rE CfejAfe (A k c k 



J\ (x) 



dx. 
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Appendix IV 
Proof of Theorem IV 

Proof: The proof of the result for the representation of the distribution and density of the 
total interference follows by taking the spherically symmetric bivariate stable distribution derived 
in Theorem [T] and utilising the closure under convolution result presented in Lemma |2] for each 
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Poisson mixture component to obtain the result in (|48)> . 

Mu) (y) = ^EP(^ = *)/Y*(«) (y) 
fc=i 

= d *L E / p = fc i A ^) ^ a ( A ^ a ' 6 ) /y^ (u) (y) d\ K 

k=l J 

= E / exp (~ Aa ") i^fTa) exp (~ bXK "> A A 7ldA ^/ Y ^ ( u) (y) dA ^ 



(84) 



*K E 

" max / j 



a + k — 1)! / o \ / 1 



£ (a-l)lfc! V 1 + fo / V! + & 

x <S Q ( sgn (c fe ) (u) , |c fc | 7 fc (u) , c fc 5 fe (u) ; 



Note, the mixture weight of the doubly stochastic compound processes (Cox process) for 
the Poisson-Gamma distribution is derived by considering Bayes Theorem and the conjugacy 
property of the Poisson-Gamma model [9] from which we know p(X\k) = T(X; a + k, f3 + l), hence: 

F(k\X)p(X) Po(k;X)T(X; aj P) 



P ^ p{X\k) T(X;a + k,f3 + l) 

^ e ~ P a ya-l ^ ( (P + -(/3+l)X 

k\ T(a) J ■ \ T(a + k) 

T{a + k) a {a + k-l)\ f f3 \ a f 1 



(85) 

k 



r(a)k\ (l + /3) Q+fc (a-l)!fc! \1 + PJ V 1 + P , 

Note, this mixing weight is then a negative binomial probability with r = a and p = as 
presented in [9]. 

Next, the representation of the density in (|50)) and distribution function in (|5T]> are obtained 
using the projection and SMiN representations in Lemma [TT] to obtain a parametric closed form 
representation which can be evaluated. ■ 
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